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^ I The process of chaiTn quark fragmentation is studied using D*^ meson production in deep- 

Q ■ inelastic scattering as measured by the HI detector at HERA. The parameters of fragmen- 

OO , tation functions are extracted for QCD models based on leading order matrix elements and 

DGLAP or CCFM evolution of partons together with string fragmentation and particle de- 
cays. Additionally, they are determined for a next-to-leading order QCD calculation in the 
^ . fixed flavour number scheme using the independent fragmentation of charm quarks to D*^ 

^ \ mesons. Two different regions of phase space are investigated defined by the presence or 

absence of a jet containing the D*^ meson in the event. The fragmentation parameters 
extracted for the two phase space regions are found to be different. 
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1 Introduction 



The production of charm quarks is expected to be well described by perturbative Quantum 
Chromodynamics (pQCD) due to the presence of a hard scale provided by the charm mass. The 
evolution of an "off-shell" charm quark via gluon radiation until it is "on-shell" can be calcu- 
lated in pQCD in fixed order of the strong coupling or by summing all orders in the leading-log 
approximation. The transition of an on-shell charm quark into a charmed hadron is, however, 
not calculable within the framework of pQCD and is thus usually described by phenomeno- 
logical models. One of the major characteristics of this transition process is the longitudi- 
nal momentum fraction transferred from the quark to the hadron, the distribution of which is 
parametrised by a fragmentation function. 

Several phenomenological models are available to describe the transition of a quark into 
hadrons, for example the independent fragmentation [1], the string [2], and the cluster model [3]. 
The fragmentation function is unambiguously defined in a given context of a phenomenological 
model together with a pQCD calculation. Universality is then only expected to hold within this 
context. 

The fragmentation function is not a directly measurable quantity as the momentum of the 
heavy quark is experimentally not accessible. Also the momentum distribution of the heavy 
hadron can only be measured within a restricted phase space. The momentum spectrum is 
further affected by the fact that some heavy hadrons are not produced directly, but are the result 
of decays of still heavier excited states, whose contribution is not well known. 

The production of charmed hadrons has been measured and parameters of fragmentation 
functions have been determined in e+e^ annihilation experiments [4]. The HI and ZEUS 
collaborations have published total cross sections for the production of various charmed hadrons 
in deep-inelastic ep scattering (DIS) [5] and in photoproduction [6]. These data show that the 
probabilities of charm quarks to fragment into various final state hadrons are consistent, within 
experimental uncertainties, for e"^e~ and ep collisions. 

In this paper the transition of a charm quark into a D*^ meson in DIS is further investigated. 
The normalised differential cross sections are measured as a function of two observables with 
different sensitivity to gluon emissions. The momentum of the charm quark is approximated 
either by the momentum of the jet, which includes the D*^ meson, or by the sum of the mo- 
menta of particles belonging to a suitably defined hemisphere containing the D*^ meson. The 
measurements are performed for two different event samples. The DIS phase space and the 
kinematic requirements on the D*^ meson are the same for both samples. In the first sample, 
referred to as the jet sample", the presence of a jet containing the D*^ meson and ex- 
ceeding a minimal transverse momentum is required as a hard scale. In the second sample no 
such jet is allowed to be present. This sample, referred to as the "no D*^ jet sample", allows 
the investigation of charm fragmentation in a region close to the kinematic threshold of charm 
production. 

The normalised differential cross sections are used to fit parameters of different fragmen- 
tation functions: for QCD models as implemented in the Monte Carlo (MC) programs RAP- 
GAP [7] and CASCADE [8], which use the Lund string model for fragmentation as imple- 
mented in PYTHIA [9], and for a next-to-leading order (NLO) QCD calculation as imple- 
mented in HVQDIS [10] with the addition of independent fragmentation of charm quarks to 
D*^ mesons. 
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The paper is organised as follows. Section [2] gives a brief description of the HI detector. 
It is followed by the details of the event selection, the D*^ meson signal extraction and the jet 
selection in section[3l The experimental fragmentation observables are defined in section|4l The 
QCD models and calculations used for data corrections and for the extraction of fragmentation 
functions are described in section [5l The data correction procedure and the determination of 
systematic uncertainties is explained in section [6l In section |7] the results of the measurements 
and of the fits of the fragmentation parameters are given. 

2 HI Detector 

The data were collected with the HI detector at HERA in the years 1999 and 2000. During 
this period HERA collided positrons of energy = 27.5 GeV with protons of energy Ep = 
920 GeV, corresponding to a centre-of-mass energy of ^/s = 319 GeV. The data sample used 
for this analysis corresponds to an integrated luminosity of 47 pb~^. 

A right handed Cartesian coordinate system is used with the origin at the nominal primary 
ep interaction vertex. The direction of the proton beam defines the positive z-axis (forward 
direction). Transverse momenta are measured in the x-y plane. Polar (9) and azimuthal (0) 
angles are measured with respect to this reference system. The pseudorapidity is defined as 
rj = — ln(tan |). 

A detailed description of the HI detector can be found in [11]. Here only the components 
relevant for this analysis are described. The scattered positron is identified and measured in 
the SpaCal [12], a lead-scintillating fibre calorimeter situated in the backward region of the HI 
detector, covering the pseudorapidity range —4.0 < rj < —1.4. The SpaCal also provides infor- 
mation to trigger on the scattered positron in the kinematic region of this analysis. Hits in the 
backward drift chamber (BDC) are used to improve the identification of the scattered positron 
and the measurement of its angle [13]. Charged particles emerging from the interaction re- 
gion are measured by the Central Silicon Track detector (CST) [14] and the Central Tracking 
Detector (CTD), which covers a range —1.74 < r] < 1.7 A. The CTD comprises two large cylin- 
drical Central Jet drift Chambers (CJCs) and two z-chambers situated concentrically around the 
beam-line, operated within a solenoidal magnetic field of 1.16 T. The CTD also provides trig- 
gering information based on track segments measured in the r-0-plane of the CJCs and on the 
2-position of the event vertex obtained from the double layers of two Multi-Wire Proportional 
Chambers (MWPCs). The tracking detectors are surrounded by a finely segmented Liquid Ar- 
gon calorimeter (LAr) [15]. It consists of an electromagnetic section with lead absorbers and a 
hadronic section with steel absorbers and covers the range — 1.5 < 77 < 3.4. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
ep — ep7, where the photon is detected in a calorimeter close to the beam pipe aX z = — 103 m. 

3 Data Selection and Analysis 

The events selected in this analysis are required to contain a scattered positron in the SpaCal and 
at least one D*^ meson candidate. The scattered positron must have an energy above 8 GeV. 
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The virtuality of the photon and the inelasticity y are determined from the measured energy 
E'^ and the polar angle 0'^ of the scattered positron via the relations: 

= AE^K cos' (I) and 2/ = 1 " f sm^ (|) • (D 
In addition, the energy W of the 7*p rest-frame is determined using: 

W'=ys-Q\ (2) 

where s = AE^Ep is the centre-of-mass energy squared of the ep system. The photon virtuality 
is required to be in the range 2 < < 100 GeV^. This kinematic range is determined 
by the geometric acceptance of the SpaCal. The inelasticity is required to lie in the region 
0.05 < y < 0.7. The difference between the total energy E and the longitudinal component Pz 
of the total momentum, as calculated from the scattered positron and the hadronic final state, is 
restricted to 40 < E — < 75 GeV. This requirement suppresses photoproduction background, 
where a hadron is misidentified as the scattered positron. It also reduces the contribution of DIS 
events with hard initial state photon radiation, where the positron or photon escapes in the 
negative ^-direction. This leads to values of E — P^ significantly lower than the expectation 
2Ee = 55 GeV. 

The D*"^ mesons are reconstructed from tracks using the decay channel D*^ D^T^t ~^ 
(jr~7r+)7r+ and its charge conjugate, where tx^ denotes the low momentum pion from the D*^ 
meson decay. Requirements on the transverse momentum and pseudorapidity of the D*^ meson 
candidate and its decay products, as well as on particle identification using AE/diX, are similar 
to those used in previous HI analyses [16]. A summary of the most important requirements is 
given in tabled 





Pt{K) > 0.25 GeV 
Prin) > 0.25 GeV 
Pt{K) + PT(vr) > 2 GeV 
\M{K7r) - M{D^) 1 < 0.07 GeV 




PT(vrs) > 0.12 GeV 

\v{D*^)\ < 1.5 

1.5 < PtP*^) < 15 GeV 



Table 1: Kinematic requirements for the selection of meson candidates. 

To select D*"^ meson candidates the invariant mass difference method [17] is used. The 
distribution of AM^.i = M{K7ms) — M{Ktx) is shown in figure [D for the full data sample, 
together with the wrong charge K^Tx^Tif combinations, using K^tx"^ pairs in the accepted 
mass range. Detailed studies show that the wrong charge /S.Md*± distribution provides a good 
description of the right charge K^-K^nf combinatorial background. 

The signal is extracted using a simultaneous fit to the ^Md*± distribution of the right and 
wrong charge combinations. The signal is fitted using a modified Gaussian function [18] 
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Gmod oc Nd*± exp [-0.5xi+i/(^+°-5^)] , (3) 

where x = | AM^i.i — MqI/o-. The signal position Mq and width a as well as the number of 
D*^ mesons Nf)*± are free parameters of the fit. The background is parametrised as a power 
function of the form N (a+1) (AMo*± —niT^Y/ (A/max— "^tt)''"*^^, with the fit boundaries given by 
the charged pion mass and M^ax = 0.17 GeV. The two free parameters a and N determine 
the shape and normalisation of the background, respectively. The total event sample is fitted 
using six free parameters: three for the modified Gaussian, two for the normalisation of the 
right and wrong charge AA//).± background distributions and one for the background shape, 
common for the right and wrong charge combinatorial background. In total 2865 ± 89 (stat.) 
D*^ mesons are obtained. For the differential distributions, the number of D*^ mesons in each 
measurement bin is extracted using the same procedure, except that the position of the signal 
peak and its width are fixed to the values determined from the fit to the total sample. 

The hadronic final state is reconstructed in each event using an energy flow algorithm. The 
algorithm combines charged particle tracks and calorimetric energy clusters, taking into account 
their respective resolution and geometric overlap, into so called hadronic objects while avoiding 
double counting of energy [19]. The hadronic objects corresponding to the three decay tracks 
forming the D*^ meson are removed from the event and replaced by one hadronic object having 
the four-momentum vector of the reconstructed D*^ meson candidate. The energy of the D*^ 
meson is calculated using M{D*^) = 2.010 GeV [20]. 

Jets are found in the 7*p rest-frame using the inclusive k-r cluster algorithm [21] with the 
distance parameter i? = 1 in the rj-cp plane. In order to combine hadronic objects into jets, the 
E-recombination scheme is applied using the four-momenta of the objects. The jet containing 
the D*^ meson candidate is referred to as the D*"^ jet and is required to have a jet transverse 
energy > 3 GeV in the 7*p rest-frameQ. According to MC simulations, the -D*^ jet is found 
to be well correlated with the original direction of the charm or anti-charm quark. The distance 
in azimuth-pseudorapidity, Ar = a/ A?7^ + A0^, between the charm quark jet, found using 
final state partons ("parton level"), and the -D*^ jet, found using final state hadrons ("hadron 
level"), is below 0.3 for 90% of all events. The correlation between the D*^ jet at hadron level 
and the D*^ jet found using charged particle tracks and calorimetric clusters ("detector level") 
is even better, since most of the energy of these jets is reconstructed from tracks, which are well 
measured in the tracking system. The number of -D*^ mesons is 1508 ± 68 (stat.) in the D*"^ 
jet sample and 1363 ± 54 (stat.) in the no D*^ jet sample. 

4 Definition of Experimental Observables 

A standard method to study fragmentation is to measure the differential production cross section 
of a heavy hadron (H) as a function of a scaled momentum or energy. In e~^e~ experiments a 
customary experimental definition of the scaled energy is Zg+e- = -EnZ-Ebcam, where -Ebcam is 
the energy of the beams in the centre-of-mass system. In leading order (LO), i.e. without gluon 

^Kinematic variables with the superscript * refer to the rest-frame of the virtual photon (7*) and proton. In this 
frame the photon direction is taken as the direction of the z-axis. The four-vector of the virtual photon used in the 
boost calculation is determined from the measurement of the scattered positron. 
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emissions, the beam energy is equal to the energy of the charm or anti-charm quark, which are 
produced in a colour singlet state. The differential cross section of heavy hadron production as 
a function of i^+e- is directly related to the fragmentation function. 

In the case of interactions the situation is more complex. In DIS the dominant process for 
D*^ meson production at HERA is photon-gluon fusion '-)*g cc [16]. In this case the cc pair 
is produced in a colour octet state. The energy of the charm quark pair depends on the energy 
of the incoming photon and gluon. Hadrons produced by initial state gluon emissions and by 
fragmentation of the proton remnant are also present in the final state. 

In this analysis charm fragmentation is studied by measuring the differential cross sections 
of meson production as a function of two different observables Zhcm and zjet defined below, 
which are sensitive to the fraction of momentum inherited by the -D*^ meson from the initial 
charm quark [22] . 



The hemisphere method: z = Zhem 

In the LO photon-gluon fusion process, which dominates charm production at HERA, the charm 
and anti-charm quarks are moving in the direction of the virtual photon in the Y'P rest-frame 
of reference. This is due to the fact that the photon is more energetic than the gluon, which 
typically carries only a small fraction of the proton's momentum. Assuming no further gluon 
radiation in the initial and final state, the charm and anti-charm quarks are balanced in transverse 
momentum (figure[2l left). This observation suggests to divide the event into hemispheres, one 
containing the fragmentation products of the charm quark, the other one those of the anti-charm 
quark. In order to suppress contributions from initial state radiation and the proton remnant, 
particles pointing in the proton direction of the 7*j> rest-frame (77* < 0) are discarded. The 
projections of the momenta of the remaining particles onto a plane perpendicular to the 7*p- 
axis are determined. Using the projected momenta, the thrust-axis in this plane, i.e. the axis 
maximising the sum of the momenta projections onto it, is found. A plane perpendicular to 
the thrust-axis allows the division of the projected event into two hemispheres, one of them 
containing the D*"^ meson and usually other particles (figure[2l right). The particles belonging 
to the same hemisphere as the D*"^ meson are attributed to the fragmentation of the charm or 
anti-charm quark. The fragmentation observable is defined as: 

where in the denominator the energy E* and the momentum P* of all particles of the D*^ 
meson hemisphere are summed. The longitudinal momentum -Pld*± defined with respect 
to the direction of the three-momentum of the hemisphere, defined as the vectorial sum of 
the three-momenta of all particles belonging to the hemisphere. The variable Zhem is invariant 
with respect to boosts along the direction of the sum of the momenta of all particles in the 
hemisphere. Neglecting the mass of the -D*^ meson and of the hemisphere, this definition of 
Zhem simplifies to the ratio of their energies. 
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The jet method: z = zjet 



In the case of the jet method the energy and direction of the charm quark are approximated by 
the energy and direction of the reconstructed jet, which contains the D*"^ meson. The fragmen- 
tation observable is defined in analogy to Zhcm as: 

where the longitudinal momentum Plj^*± is defined with respect to the direction of the three- 
momentum of the jet. The jet finding and the determination of zjct are performed in the YP 
rest-frame. 

Both fragmentation observables are defined in such a way that they would lead to similar 
distributions, assuming independent fragmentation and no gluon radiation. The measured dis- 
tributions, however, are expected to differ, as they have different sensitivities to gluon radiation 
and charm quarks, which are colour connected to the partons of the proton remnant. The hemi- 
sphere method typically includes more energy around the charm quark direction than the jet 
method. The parameters of fragmentation functions should however be the same, if extracted 
for a QCD model, which provides a very good description of the underlying physics over the 
full phase space of this analysis. A comparison of both methods thus may provide a consistency 
check and a test of the perturbative and non-perturbative physics as encoded in the models. 

The measurement is restricted to the regions 0.2 < Zhcm < 10 and 0.3 < zjct < 10, 
as at lower z it is not possible to separate the D*^ meson signal from the large combinatorial 
background. In order to minimise the sensitivity of the analysis to the total -D*^ meson cross 
section, and to reduce systematic errors, normalised differential cross sections are measured as 
a function of the fragmentation observables Zhem and zjet- The normalisations are chosen such 
that their integrals over the respective z-regions yield unity. 



5 QCD Models and Calculations 

The MC programs RAPGAP and CASCADE are used to generate events containing charm and 
beauty quarks, which are passed through a detailed simulation of the detector response, based on 
the GEANT simulation program [23]. They are reconstructed using the same software as used 
for the data. These event samples are used to determine the acceptance and efficiency of the 
detector and to estimate the systematic errors associated with the measurements. In addition, 
these models are fitted to the data in order to determine the parameters of the fragmentation 
functions. 

The Monte Carlo program RAPGAP [7], based on coUinear factorisation and DGLAP [24] 
evolution, is used to generate the direct process of photon-gluon fusion to a heavy cc pair, 
where the photon acts as a point-like object. In addition, RAPGAP allows the simulation of 
charm production via resolved processes, where the photon fluctuates into partons, one of which 
interacts with a parton in the proton, and the remaining partons produce the photon remnant. 
The program uses LO matrix elements with massive (massless) charm quarks for the direct 
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(resolved) processes. Parton showers based on DGLAP evolution are used to model higher 
order QCD effects. 

The CASCADE program [8] is based on the /cT-factorisation approach. Here, the calcula- 
tion of the photon-gluon fusion matrix element takes into account the charm quark mass and 
the virtuality and the transverse momentum of the incoming gluon. Gluon radiation off the 
incoming gluon as well as parton showers off the charm or anti-charm quark are implemented 
including angular ordering constraints. The gluon density of the proton is evolved according to 
the CCFM equations [25]. The fcx-unintegrated gluon density function AO [26], extracted from 
inclusive DIS data, is used. 

In both RAPGAP and CASCADE the hadronisation of partons is performed using the Lund 
string model as implemented in PYTHIA [9]. In the Lund model, the heavy hadron is produced 
in the process of string breaking. The fraction of the string longitudinal momentum z carried 
by the hadron is generated according to different choices of adjustable fragmentation functions 
Dq(2;). Within this analysis three widely used parametrisations are employed, of which two 
depend on a single free parameter, and one depends on two free parameters. The parametrisation 
suggested by Peterson et al. [27] has the functional form: 



D'^[z\ oc — — ——^ — — , (6) 

^ z[l - {l/z)-e/{l - z)Y ' ^ ^ 



and the one by Kartvelishvili et al. [28] is given by: 



D^{z)<xz"{l-z) . (7) 

The free parameters e and a determine the "hardness" of the fragmentation function and are 
specific to the flavour of the heavy quark, i.e. charm in the case of D*"^ meson production. The 
parametrisation inspired by Bowler and Morris [29] (referred to as the Bowler parametrisation) 
has the functional form: 



^Sw«7;;^(i--r-p(-^)- (8) 

The shape of the fragmentation function is determined by two free parameters a and b, uiq is 
the mass of the heavy quark, Mt = ^/ + P| the transverse mass of the heavy hadron, and 
rq = 1 as defauh in PYTHIA. 

For data corrections the parameter setting tuned by the ALEPH collaboration [30] together 
with the Peterson fragmentation function is used for the fragmentation of partons in PYTHIA. 
It includes higher excited charm states, of which some also decay to D*"^ mesons and contribute 
significantly to the D*"^ meson yield. When extracting parameters of the fragmentation func- 
tions also the default parameter setting of PYTHIA is used as an alternative. In this case no 
higher excited charm states are produced. Both parameter settings are indicated in table [21 

The parameters for the Kartvelishvili and Peterson fragmentation functions are also ex- 
tracted for the HVQDIS program [10]. HVQDIS is based on the NLO, i.e. 0{a'^), calculation 
in the fixed flavour number scheme, with three light active flavours as well as gluons in the pro- 
ton. The proton parton density functions (PDFs) of the light quarks and the gluon are evolved 
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according to the DGLAP equations. Massive charm quarks are assumed to be produced only 
perturbatively via photon-gluon fusion and higher order processes. The final state charm quarks 
are fragmented independently into D*^ mesons in the rest-frame. Kartvelishvili and Pe- 
terson parametrisations are used to generate the charm quark's momentum fraction transferred 
to the D*"^ meson. The energy of the charm quark is calculated using the on-mass-shell con- 
dition. In addition, the D*"^ meson can be given a transverse momentum Pt with respect to 
the charm quark, according to the function Px exp(— /3Pt). The value used for the parameter P 
corresponds to an average Pt(-D*^) of 350 MeV. 

The Monte Carlo programs RAPGAP and HERWIG [31] are used to estimate the size of 
the hadronisation corrections to the data for comparison with HVQDIS predictions. While the 
perturbative QCD model of HERWIG is similar to the one of RAPGAP, the HERWIG program 
employs the cluster hadronisation model, which is quite different from the Lund string model 
used by PYTHIA. 

The basic parameter choices for the QCD models and the NLO calculation are summarised 
in tabled 



6 Data Corrections and Systematic Errors 

In this analysis, the differential cross section for the production of D*^ mesons, which result 
from the fragmentation of charm quarks either directly or via decays from higher excited charm 
states, is measured. The small contribution of P**^ mesons originating from B-hadron decays is 
estimated with RAPGAP and is subtracted from the data. It is on the level of 1 to 2%. The data 
are corrected for detector and QED radiative effects. The transverse momentum and pseudo- 
rapidity distributions of the P*^ mesons of the Monte Carlo event samples, which are used to 
correct the data samples for detector effects, are reweighted to the corresponding distributions 
of the data in order to achieve an improved description. The rj and Pt dependent reweighting 
factors differ from unity by typically 10 — 30%. After this reweighting, both RAPGAP and 
CASCADE provide a good description of the data as shown in figure [H The description of the 
no P*^ event sample by the reweighted MC models, as shown in figure IH is worse. 

The measurement bins are defined in such a way that the purity in each bin, defined as the 
fraction of events reconstructed in a Zhcm or zjet bin that originate from that bin on hadron level, 
is between 40 and 70%. 

The correction for the detector effects is done using regularised deconvolution, taking into 
account migrations between measurement bins [32]. The detector response matrix is generated 
using RAPGAP, and the value of the regularisation parameter is determined through decompo- 
sition of the data into eigenvectors of the detector response matrix. As a check, the detector 
response matrix was also generated using CASCADE and found to be consistent with the one 
from RAPGAP. Statistical errors are calculated by error propagation using the covariance ma- 
trix. The data are then corrected for migrations into the visible phase space using RAPGAP 
and CASCADE. The effects of QED radiation are corrected for using the HERACLES [33] 
program, which is interfaced with RAPGAP. Correction factors are calculated from the ratio 
between cross sections obtained from the model including and not including QED radiation. 



11 



The corrections are applied bin-by-bin in Zhom and zjet- In the case of zjct the corrections are 
about 2%. In the case of Zhem^ where photons radiated into the D*^ meson hemisphere can 
be mistaken as fragmentation products of the charm quark, the corrections reach 10% for the 
lowest value of Zhcm- 

In contrast to the QCD models discussed so far, the HVQDIS program provides only a par- 
tonic final state with the exception of the additional D*"^ meson from charm fragmentation. In 
order to compare the NLO predictions to the measurements, the data are corrected to the parton 
level by means of hadronisation corrections, which are estimated using the MC generators RAP- 
GAP and HERWIG. While the quantity {E* + P£)d.± in equation Hand [5] is calculated using 
the momentum of the D*"^ meson, the jet finding and the calculation of the jet and hemisphere 
quantities, the denominators in equations |4] and [51 are performed using the partonic final state. 
All partons after parton showering are considered, and the same jet and hemisphere finding 
algorithms are applied at parton and hadron level. For each z-bin the hadronisation correction 
factor is calculated as the ratio of parton to hadron level cross section. The arithmetic mean of 
the hadronisation correction factors of both models is used to multiply the data cross section. 
In case of Zhcm the hadronisation corrections differ from unity by typically ±40%. For zjct they 
differ from unity by typically ±20%, except for the highest z-bin, where they are about 50%. 
The hadronisation corrections as determined by HERWIG and RAPGAP are similar for most 
of the measurement bins, with the exception of the lowest bin in zjet, where they differ by about 
60%. 

The following systematic uncertainties on the normalised differential cross sections are con- 
sidered: 

• The energy uncertainty of the scattered positron varies linearly from ±3% for an energy 
of 8 GeV to ±1% for 27 GeV. 

• The polar angle of the scattered positron has an estimated uncertainty of ±1 mrad. 

• The uncertainty of the energy scale of the hadronic objects is made up of ±0.5% due to 
tracks and ±4% (±7%) due to LAr (Spacal) clusters. 

• The effect of the uncertainty of the tracking efficiency on reconstructing the D*"^ meson is 
determined by changing the nominal efficiency in the simulation as a function of track i] 
and Pt- In the central region of the accepted tj-Pt phase space the estimated uncertainty 
of the nominal efficiency is ±2%, in the regions of large |?7| but not small Pt it is ±3%, 
and for large |?7| and small Pt it is ±4%. 

• The value of dE/dx of the -D*^ meson decay products has an estimated uncertainty of 
±8%, which is of similar size as the experimental resolution in dE/dx. 

• The uncertainty of the D*"^ meson signal extraction is estimated using different D*"^ 
counting techniques and by using different fit functions for the background parametrisa- 
tion. The largest uncertainty comes from the background description, which determines 
the systematic error on the signal extraction. 

• The uncertainty of beauty production by the RAPGAP MC is assumed to be ±100%. The 
resulting small uncertainty of the normalised D*"^ meson cross sections is taken to be 
symmetrical. 
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• The effect of using different MC models for the small correction for migrations into the 
visible phase space is studied using RAPGAP and CASCADE. The factors used to correct 
the data are determined as the average of the correction factors obtained from the two 
models. Half the difference is taken as systematic uncertainty. 

• For parton level corrected distributions half the difference between the hadronisation cor- 
rection factors of RAPGAP and HERWIG is taken as the uncertainty due to the different 
fragmentation models. 

Other systematic effects, which are investigated and found to be negligible, are: the effect 
of reflections, i.e. wrongly or incompletely reconstructed D*"^ meson decays, on the shape 
of the fragmentation observables, the effect on acceptance and reconstruction efficiency from 
including diffractive events, the effect of using different MC models for the deconvolution of 
the data and the uncertainty of the QED radiative effects. 

Each source of systematic error is varied in the Monte Carlo simulation within its uncer- 
tainty. In each measurement bin, the corresponding deviation of the normalised cross sections 
from the central value is taken as the systematic error. Among the systematic errors the un- 
certainties due to the scattered positron energy scale, the hadronic energy scale, and the beauty 
fraction are correlated amongst the bins in z. In the extraction of the parameters of the frag- 
mentation functions, the statistical and systematic errors with their correlations are taken into 
account. The average effect of various systematic errors on the Zhem and zjot distributions is 
summarised in table IH Since the distributions of Zhom and zjet are normalised, the effect of 
many systematic uncertainties is reduced and the statistical error dominates the uncertainty of 
the measurement. 



7 Results 

7.1 Normalised differential cross sections and comparison with predic- 
tions 

The differential cross sections of -D*^ meson production as a function of the fragmentation 
observables Zhcm and zjct are shown in figure[5]for the D*"^ jet sample. They refer to the visible 
phase space given by 2 < < iqo GeV^, 0.05 < y < 0.7, 1.5 < Pt{D*^) < 15 GeV and 
|?7(-D*^)| < 1.5. In addition, a D*"^ jet with > 3 GeV in the YP rest-frame is required 
in order to have the same hard scale in the event for both distributions Zhcm and zjct- The 
measurements and the corresponding predictions are normalised such that their integrals over 
the respective z-regions yield unity. The normalised cross sections and their errors are given in 
table [5] for the hemisphere observable and in table [6] for the jet observable. 

Figure [5] also includes predictions of RAPGAP with three commonly used fragmentation 
parameter settings for PYTHIA (described in table O, obtained from e'^e~ annihilation. The 
settings and the corresponding values of x^/n.d.f., as calculated from the data and the model 
predictions, are summarised in table |7l In general, there is reasonable agreement between data 
and the QCD model with all settings for both the jet and the hemisphere observables. The 
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large difference between the two distributions observed in the highest zjet bin is mainly due to 
a significant fraction of D*^ jets consisting of a D*^ meson only, for which zjot equals unity. 
CASCADE provides a similar description of the data as RAPGAP. 

7.2 Extraction of parameters for the Kartvelishvili and Peterson frag- 
mentation functions 

The normalised D*^ meson differential cross sections as a function of Zhcm and zjct are used 
to extract optimal parameters for the Peterson and Kartvelishvili fragmentation functions de- 
scribed in section [51 Both parametrisations have a single free parameter. 

The parameter extraction is done by comparing different model configurations with the data. 
A configuration is defined by one of the QCD calculations (RAPGAP, CASCADE or HVQDIS), 
by one of the fragmentation functions (Peterson or Kartvelishvili) and by a possible value for the 
corresponding fragmentation parameter, e ora. For RAPGAP and CASCADE the configuration 
also depends on the PYTHIA parameter settings used (ALEPH and default, see tableO. In order 
to be able to compare all configurations to the data, a reweighting procedure is applied. For each 
of the QCD calculations large event samples with D*"^ mesons are generated using the Peterson 
fragmentation function. For these events the 2;-value of the fragmentation function, used by 
the model to generate the fraction of charm quark or string momentum transferred to the D*"^ 
meson, is stored such that each event can be reweighted to another fragmentation function or 
any other parameter value. For each configuration the predicted and measured distributions of 
the fragmentation observables are used to determine a as a function of the fragmentation 
parameter. In the calculation of the the full covariance matrix is used, taking into account 
correlated and uncorrelated statistical and systematic errors. The fragmentation parameter is 
determined at the minimum of the x^- The shape of the distribution is used to determine the 
±lcr error (using Xmm + 1) of the extracted parameter. As an example, in figure [6] the data are 
compared to the prediction of RAPGAP with the ALEPH setting for PYTHIA as given in table[2l 
but using the Kartvelishvili parametrisation. The two lines indicate the ±la total uncertainty 
around the best fit value of a. The description of the data by CASCADE is similar. 

The parameters a and e, which are extracted using RAPGAP and CASCADE, with and 
without higher excited charmed hadrons, are summarised in table [8] together with their corre- 
sponding values of x^/n.d.f.. With the fitted parameters the model predictions using either the 
Peterson or the Kartvelishvili parametrisations describe the data reasonably well, with the Kar- 
tvelishvili parametrisation being in all cases slightly preferable, as indicated by the values of 
X^/n.d.f.. When using the same PYTHIA parameter setting, the fragmentation parameters a 
and e, extracted from the Zhcm and Zjet observables, are in good agreement. Both RAPGAP and 
CASCADE lead to statistically compatible parameters e and a. A priori, agreement in the frag- 
mentation function parameters for RAPGAP and CASCADE is not required, since the models 
differ in terms of simulated processes (direct and resolved in case of RAPGAP compared to 
direct only for CASCADE) and in their implementation of perturbative QCD. 

The fragmentation parameters a and e depend significantly on the PYTHIA parameter set- 
tings used, i.e. whether D*^ mesons are assumed to be produced only via direct fragmentation 
of charm quarks or additionally originate from decays of higher excited charm states. In the 
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latter case the D*"^ mesons carry a smaller fraction of the original charm or anti-charm quark 
momentum in comparison with the directly produced ones. Both the default PYTHIA setting 
and the setting containing higher excited charm states describe the data equally well. The values 
of the Peterson parameter e extracted for the PYTHIA setting containing higher charm states, 
see table [HI are in agreement with the value e = 0.04 tuned by ALEPH [30]. This result is 
consistent with the hypothesis of fragmentation universality in ep and e+e~ processes. 

The NLO calculation as implemented in HVQDIS with the Kartvelishvili fragmentation 
function leads to a good fit of the data, corrected for hadronisation effects, as shown in figure |7l 
On the other hand HVQDIS provides a rather poor description of Zhem and zjct, if the Peterson 
fragmentation function is used (the values of the fit are shown in table [8]). Simulating a Pt 
of the D*"^ meson with respect to the charm quark direction, as explained in section[5l has only 
a little effect on the extracted value of a. 

The distributions of Zhcm and zjct are also measured in two bins of and W. In the case 
of the bins are defined as 2 < < ig QeV^ and 10 < < iqo GeV^. The accessible 
range in W, determined by the cuts on and y, corresponds to 70 < 11^ < 270 GeV, and the 
bins are defined as 70 < < 170 GeV and 170 < < 270 GeV. Correction factors and 
systematic uncertainties for these samples are determined in the same way as for the full D*"^ 
jet sample. The data are compared to the QCD models with the PYTHIA parameter setting 
including higher excited charm states and using the Kartvelishvili fragmentation function. For 
the low and high bins the measured distributions are found to be almost the same and well 
described by the QCD models. The distributions of zjct for the low and high W regions are also 
similar. A difference is observed for the Zhcm distribution, which is softer at high W as shown 
in figure [8l RAPGAP and CASCADE show the same behaviour as a function of W as observed 
in data. This behaviour can be understood as being partly due to enhanced gluon radiation at 
high W and partly due to the kinematic effect of the requirement Pt(-D*^) > 1.5 GeV. For 
events at low W, where the charm quark tends to have smaller energy than at high W, a D*"^ 
meson needs to carry a large fraction of the original quark momentum in order to pass the Pt 
requirement. 

The hemisphere observable allows an investigation of charm fragmentation close to the 
kinematic threshold, at the limit of applicability of the concept of fragmentation functions, by 
selecting events without a -D*^ jet with > 3 GeV. As estimated by MC, the mean centre- 
of-mass energy squared of the 'j*g system s for this sample is about 36 GeV^, to be compared 
with about 100 GeV^ for the D*"^ jet sample. This event sample (no -D*^ jet sample) has no 
overlap with the D*^ jet sample investigated so far. The normalised cross section as a function 
of Zhcm for the no D*^ jet sample is shown in figure |9] and listed in table |9l Predictions of 
RAPGAP with the three commonly used fragmentation parameter settings for PYTHIA (see 
table [21), which provide a reasonable description of the D*"^ jet sample (see table [7]), fail for 
this sample. Also the prediction using the fragmentation parameters obtained from the D*"^ jet 
sample is not able to describe these data. 

The fragmentation parameters for RAPGAP, CASCADE and the NLO calculation are ex- 
tracted from the no D*^ jet sample using the same procedure as for the D*"^ jet sample. The fit 
results are summarised in table [TOj The predictions of RAPGAP, showing the ±la total uncer- 
tainty around the fitted value of a are also presented in figure [9l The fragmentation parameters 
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obtained for RAPGAP and CASCADE are statistically compatible. The fragmentation param- 
eters fitted to the no D*^ jet sample are found to be significantly different from those for the 
D*^ jet sample. They indicate that the fragmentation function for an optimal description of 
the sample without a D*"^ jet needs to be significantly harder than for the D*"^ jet sample. The 
NLO calculation as implemented in HVQDIS fails to describe the no D*^ jet sample as shown 
in figure [TOl 

Several parameters of the QCD models, for example those influencing parton showers, have 
been varied, in trying to describe both samples using the same value for the fragmentation func- 
tion parameter. However, it was not possible to find MC parameters leading to a consistent 
fragmentation function for the two samples. Furthermore, the effect of diffractive production 
of D*"^ mesons was not able to explain the difference between the fragmentation parameters 
observed for the two samples. These investigations indicate that QCD models, together with 
simple parametrisations of the fragmentation functions, are not able to describe charm frag- 
mentation consistently in the full phase space down to the kinematic threshold. 



8 Conclusions 

The fragmentation of charm quarks into D*"^ mesons in DIS is studied using the HI detector 
at the HERA collider. The normalised D*^ meson differential cross sections are measured as 
a function of two observables sensitive to fragmentation, the hemisphere observable Zhem and 
the jet observable zjct, in the visible DIS phase space defined by 2 < < 100 GeV^ and 
0.05 < y < 0.7, and the D*^ meson phase space defined by 1.5 < P-r{D*^) < 15 GeV and 
|?7(D*^)| < 1.5. An additional jet with > 3 GeV, containing the meson, is required in 
the 7*p rest-frame in order to provide a hard scale for the events. 

The data are compared with predictions of RAPGAP with three widely used PYTHIA pa- 
rameter settings and the Peterson and the Bowler parametrisations for the fragmentation of 
heavy flavours obtained from e^e^ annihilation. They provide a reasonable description of the 
ep data presented. 

The normalised differential cross sections are used to fit the parameters of the Kartvelishvili 
and Peterson fragmentation functions within the framework of the QCD models RAPGAP and 
CASCADE. The fragmentation parameters extracted using the Zhem and zjct observables are in 
good agreement with each other. Both QCD models lead to statistically compatible parameters. 
The value of the Peterson parameter e extracted for the PYTHIA parameter setting, which 
includes not only D*"^ mesons from direct fragmentation of charm quarks but also from the 
decays of higher excited charm states, is in agreement with the value of e = 0.04 tuned by 
ALEPH. This result is consistent with the hypothesis of fragmentation universality between ep 
and e~^e~ collisions. 

The QCD models, with the fragmentation parameters fitted to the data, also provide a good 
description of the and W dependence of the fragmentation observables. 

The data, corrected to the parton level, are also compared to the NLO calculation as imple- 
mented in HVQDIS, with the addition of independent fragmentation of charm quarks to D*"^ 
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mesons. A good fit to the data is obtained when using the fragmentation function by Kartveli- 
shvili et al., while using the one of Peterson et al. results in a poor fit. 

Finally, the hemisphere method is used to study the fragmentation of charm produced close 
to the kinematic threshold, by selecting a data sample fulfilling the nominal requirements on 
the DIS and D*^ meson phase space, but without a D*"^ jet having > 3 GeV in the event. 
The fragmentation parameters extracted for the QCD models, using this sample of events, are 
significantly different from those fitted to the D*"^ jet sample. Furthermore, the fit for the 
NLO calculation using the no D*^ jet sample fails. Both observations can be interpreted as an 
inadequacy of the QCD models and the NLO calculation to provide a consistent description of 
the full phase space down to the kinematic threshold. 
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PYTHIA 


ALEPH 


Default 


Description 


parameter 


setting 


setting 




MSTJ(12) 


2 


2 


baryon model option 


MSTJ(46) 





3 


parton shower azimut. corr. 


MSTJ(51) 








Bose-Einstein correlations off 


PARJ(l) 


0.108 


0.100 


P(qq)/P(q) 


PARJ(2) 


0.286 


0.300 


P(s)/P(u) 


PARJ(3) 


0.690 


0.400 


P(us)/P(ud)/P(s)/P(d) 


PARJ(4) 


0.050 


0.050 


(l/3)P(ud_l)/P(ud_0) 


PARJ(ll) 


0.553 


0.500 


P(S=l)d,u 


PARJ(12) 


0.470 


0.600 


P(S=l)s 


PARJ(13) 


0.650 


0.750 


P(S=l)c,b 


PARJ(14) 


0.120 


0.000 


P(S=0,L=1,J=1) AXIAL 


PARJ(15) 


0.040 


0.000 


P(S=1,L=1,J=0) SCALAR 


PARJ(16) 


0.120 


0.000 


P(S=1,L=1,J=1) AXIAL 


PARJ(17) 


0.200 


0.000 


P(S=1,L=1,J=2) TENSOR 


PARJ(19) 


0.550 


1.000 


extra baryon suppression 


PARJ(21) 


0.366 


0.360 


On 


PARJ(25) 


1.000 


1.000 


extra t] suppression 


PARJ(26) 


0.276 


0.400 


extra r( suppression 


PARJ(41) 


0.400 


0.300 


Lund symm. fragm.: a 


PARJ(42) 


0.885 


0.580 


Lund symm. fragm.: b 


PARJ(54) 


-0.040 


-0.050 


Peterson fragm.: —e^ 


PARJ(55) 


-0.002 


-0.005 


Peterson fragm.: —Sh 


PARJ(82) 


1.390 


1.000 


Qo 


PARP(72) 


0.295 


0.250 


A for as in time-like parton 
showers 



Table 2: PYTHIA (version 6.2) parameter settings: ALEPH [30] and defauh. The ALEPH 
setting together with the Peterson fragmentation function is used for detector corrections. 
Detailed explanation of the parameters can be found in [9] . 
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RAPGAP 


CASCADE 


HERWIG 


HVQDIS 


Proton PDFs 


CTEQ5L [34] 


AO [26] 


CTEQ5L [34] 


CTEQ5F3 [34] 


Photon PDFs 


SaSD-2D [35] 




SaSG-lD [35] 




Renorm. scale fir 


v/g2 + Pi' 








Fact, scale /i/ 




+ Qt 




y/Q^ + Ami 


Fragmentation model 


Lund string 


Lund string 


cluster 


independent 



Table 3: Parton density functions (PDFs), fragmentation models and basic parameters used in 
the QCD models and the NLO calculation. The mass nic of the charm quark is 1.5 GeV in all 
cases. The transverse momentum of the charm quark in the 7*j9 rest-frame is given by P^. The 
invariant mass squared and the transverse momentum squared of the cc pair are denoted by s 
and Q'^, respectively. 





D*^ jet sample 


No D*^ jet sample 




Zhem error 


Zjct error 


Zhem error 


Statistical uncertainty 


9.5% 


10.9% 


10.9% 


Source of systematic uncertainty 








Scattered positron energy scale 


0.8% 


0.5% 


0.5% 


Positron scattering angle 


0.1% 


0.1% 


0.1% 


Hadronic energy scale 


3.0% 


2.5% 


2.1% 


Track reconstruction efficiency 


0.1% 


0.1% 


0.1% 


dE/diX measurement 


0.1% 


0.3% 


0.8% 


D*^ signal extraction 


3.0% 


3.0% 


2.3% 


Beauty fraction 


1.2% 


0.9% 


0.6% 


Migrations into the visible phase space 


0.1% 


0.3% 


0.1% 


Total systematic uncertainty (hadron level) 


4.6% 


4.2% 


3.5% 


Hadronisation effects 


3.9% 


9.6% 


2.4% 


Total syst. uncertainty (parton level) 


6.3% 


11.1% 


4.6% 



Table 4: Experimental and model systematic uncertainties of the normalised z distributions, 
averaged over all bins. The last two uncertainties in the table apply only when data are addition- 
ally corrected for hadronisation effects to be compared with HVQDIS. The table also provides 
the statistical errors, averaged over all bins, for comparison with the systematic uncertainties. 



21 



D*^ jet sample: Zhcm 


Bin in Zhom 


1 do- 
(J dz^^(,^ 


Statistical 
eiTor 


Uncorrelated 
systematic error 


CoiTelated systematic enors 
Positron energy Hadronic scale 


Beauty 


Total eiTor 


[0.2 - 0.4[ 


0.93 


0.11 


0.03 


-0.008 


-0.027 


-0.020 


0.12 


[0.4-0.5[ 


1.53 


0.13 


0.05 


-0.014 


-0.037 


-0.019 


0.15 


[0.5 - 0.625[ 


1.80 


0.15 


0.05 


+0.002 


-0.032 


+0.007 


0.16 


[0.625 - 0.75[ 


1.85 


0.14 


0.06 


+0.008 


+0.030 


+0.019 


0.16 


[0.75 - 0.85[ 


1.27 


0.11 


0.04 


+0.008 


+0.056 


+0.016 


0.13 


[0.85 - 1.0] 


0.52 


0.06 


0.02 


+0.010 


+0.025 


+0.007 


0.07 



Table 5: Normalised _D*^ meson differential cross sections as a function of Zhcm for the D*"^ 
jet sample, in the visible phase space described in section I7.1[ The measurements are nor- 
malised such that their integral over the Zhcm range yields unity. All errors are considered to be 
symmetric in each bin. For correlated systematic errors a relative sign is indicated. 



D*^ jet sample: zjct 


Bin in zj^t 


1 do- 
o- dzjet 


Statistical 
error 


Uncorrelated 
systematic error 


CoiTelated systematic en'ors 
Positron energy Hadronic scale 


Beauty 


Total eiTor 


[0.3 - 0.55[ 


0.61 


0.10 


0.02 


-0.005 


-0.029 


-0.016 


0.11 


[0.55 - 0.7[ 


1.76 


0.15 


0.05 


-0.001 


-0.026 


+0.004 


0.16 


[0.7-0.825[ 


2.17 


0.18 


0.07 


+0.009 


+0.032 


+0.017 


0.20 


[0.825 - 0.9[ 


1.47 


0.18 


0.05 


-0.011 


+0.043 


+0.004 


0.19 


[0.9 - 1.0] 


2.03 


0.17 


0.06 


+0.012 


+0.038 


+0.011 


0.19 



Table 6: Normalised D*"^ meson differential cross sections as a function of zjct for the -D*^ jet 
sample, in the visible phase space described in section ITTl The measurements are normalised 
such that their integral over the zjct range yields unity. All errors are considered to be symmetric 
in each bin. For correlated systematic errors a relative sign is indicated. 





sam 


jet 
)le 


No D*^ jet 
sample 


RAPGAP using fi 
Parameter setting 


ragmentation by PYTHIA 
Fragmentation function 


Hemisphere 
(x'/n.d.f.) 


Jet 
(xVn.d.f.) 


Hemisphere 
(x'/n.d.f.) 


Aleph 
Default 
Default 


Peterson £ = 0.04 
Peterson e = 0.05 
Bowler a = 0.3,6 = 0.58 


5.1/5 
5.3/5 
4.8/5 


4.0/4 
5.8/4 
3.3/4 


34.2/5 
30.0/5 
21.1/5 



Table 7: The PYTHIA parameter settings and fragmentation functions used for the RAPGAP 
predictions and the corresponding values of x^/n.d.f. for the Z)*^ jet as well as the no D*"^ jet 
data samples, in the visible phase space described in section ITTl 



22 



D*"^ jet sample 


Model 


a Kartvelishvili 
(x'/n.d.f.) 


e Peterson 

(x'/n.d.f.) 




Hemisphere 


Jet 


Hemisphere 


Jet 


PYTHIA default parameter setting: 


RAPGAP 


a = 3.3+°:^ 
(1.6/4) 


« = s.i^SJ 

(2.2/3) 


(5.3/4) 


^ _ Q nfii +0.011 
(4.2/3) 


CASCADE 


a = 3.5t'oi 
(2.1/4) 


a = 3.2trs 
(3.2/3) 


^ _ n n4c;+0.012 
fc — U.U4d_o.o09 

(5.7/4) 


^ _ Q nRn+0.011 

fc — U.UUU_o.o09 

(4.7/3) 


PYTHIA with ALEPH parameter setting: 


RAPGAP 


a = 4.4i°| 
(3.0/4) 


a = 4.3l°i 
(2.8/3) 


fc — U.UOU_o.006 

(4.0/4) 


(3.8/3) 


CASCADE 


« = 4.5^°:^ 
(2.4/4) 


a = 4.4l0i 
(2.4/3) 


(3.3/4) 


^ _ n nQ4+0.007 
(3.5/3) 


Fixed-order (NLO) calculation: 


HVQDIS 


a = S.St'oi 
(4.4/4) 


a = 3.8tli 
(4.9/3) 


e = OMStlfi 
(18.3/4) 


^ _ n nQ4+0.004 
(23.3/3) 



Table 8: Fragmentation function parameters extracted for the QCD models RAPGAP and 
CASCADE, with the PYTHIA parameter settings as summarised in table [2l and for the NLO 
QCD program HVQDIS, using the hemisphere and jet observables measured with the D*"^ jet 
sample in the visible phase space described in section [TTJ 









No D*^ 


jet sample: Zhem 








Bin in Zhem 


1 dCT 


Statistical 
error 


Uncorrelated 
systematic error 


CoiTelated systematic en'ors 
Positron energy Hadronic scale 


Beauty 


Total eiTor 


[0.2 - 0.4[ 


0.50 


0.09 


0.01 




+0.003 


-0.017 


-0.007 


0.09 


[0.4-0.5[ 


0.97 


0.12 


0.02 




+0.004 


-0.023 


-0.009 


0.12 


[0.5 - 0.625[ 


1.44 


0.16 


0.03 




-0.006 


-0.026 


-0.002 


0.16 


[0.625 - 0.75[ 


1.77 


0.17 


0.04 




-0.016 


-0.005 


+0.005 


0.18 


[0.75 - 0.85[ 


2.13 


0.15 


0.05 




+0.008 


+0.037 


+0.009 


0.16 


[0.85 - 1.0] 


1.26 


0.10 


0.03 




+0.006 


+0.038 


+0.006 


0.11 



Table 9: Normalised D*"^ meson differential cross sections as a function of Zhcm for the no 
D*"^ jet sample, in the visible phase space described in section 17.11 The measurements are 
normalised such that their integral over the Zhcm range yields unity. All errors are considered to 
be symmetric in each bin. For correlated systematic errors a relative sign is indicated. 
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No D*^ jet sample 




Hemisphere observable 


Model 


a Kartvelishvili 
(x'/n.d.f.) 


e Peterson 

(x'/n.d.f.) 


PYTHIA default parameter setting: 


RAPGAP 


a = 7.5tli 
(5.5/4) 


p. -n nin+0.003 
fc — U.UiU_QQ03 

(3.9/4) 


CASCADE 


a = 6.9tli 
(4.1/4) 


^ _ 1 4+0.004 

fc — U.Ui^_QQQ3 

(2.9/4) 


PYTHIA with ALEPH parameter setting: 


RAPGAP 


a = 10.3t\i 
(2.9/4) 


fc — U.UUU_QQQ2 

(1.6/4) 


CASCADE 


a = 8At\i 
(4.6/4) 


^ _ Q nin+0.003 
(4.1/4) 


Fixed-order (NLO) calculation: 


HVQDIS 


a = 6.11°| 
(37.6/4) 


^ _ n nn7+o-ooi 

(38.6/4) 



Table 10: Fragmentation function parameters extracted for the QCD models of RAPGAP and 
CASCADE, with PYTHIA parameter settings as summarised in tabled and for the NLO QCD 
program HVQDIS, using the hemisphere observable measured with the no D*"^ jet sample in 
the visible phase space described in section [TTJ 
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Figure 1: Distributions of AA//p.± = M^Kmis) — M{Kti) for right charge combinations 
{K^Tx^Txf) and for wrong charge (i^=''7r=''7r^) combinations in the accepted mass window. 



D* hemisphere / 




Figure 2: Illustration of the hemisphere method: a cc pair and the products of its fragmentation 
in the 7*p rest-frame (left) and in a plane perpendicular to the photon momentum (right). 
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D*^ jet sample 




E\(D*jet) [GeV] ri*(D*iet) 



Figure 3: Comparison on detector level between the D* jet data sample and the reweighted 
Monte Carlo models (see section[6l) used to correct the data for detector effects. Shown are 
and 7]* of the D*"^ meson hemisphere, calculated from the sum of momenta of all particles in the 
hemisphere, and and 77* of the D*"^ jet. All observables are calculated in the rest-frame. 
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No D*^ jet sample 



> 

O 0-4 

0.3 

* 

o 

0.2 

* 

lU 

1 



m 



f H1 Data 

— RAPGAP 

- - CASCADE 



5 10 
E\(D*hem) [GeV] 




Ti*(D*hem) 



Figure 4: Comparison on detector level between the no D*^ jet sample and reweighted Monte 
Carlo models (see section (6]) used to correct the data for detector effects for the no D*^ jet 
sample. Shown are and i]* of the D*^ meson hemisphere, calculated from the sum of 
momenta of all particles in the hemisphere. All quantities are calculated in the 7*p rest-frame. 
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Figure 5: Normalised meson cross sections as a function of zjct and Zhcm for the -D*^ 
jet sample. The measurements are normalised to unity in the displayed range of zjct and Zhcm, 
respectively. The data are compared with MC predictions of RAPGAP, using PYTHIA default 
settings with Peterson or Bowler parametrisations and the ALEPH setting, which includes the 
production of higher excited charm states (see table [21). The ratio R = MC/data is shown 
as well as the relative statistical uncertainties (inner error bars) and the relative statistical and 
systematic uncertainties added in quadrature (outer error bars) for the data points put to R = 1. 
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Figure 6: Normalised D*"^ meson cross sections as a function of zjot and Zhcm for the D*^ 
jet sample. The measurements are normalised to unity in the displayed range of zjct and Zhcm, 
respectively. The same data as in figure [5] are compared to the predictions of the MC program 
RAPGAP with the ALEPH setting for PYTHIA and Kartvelishvili parametrisation using the 
fragmentation function parameter a fitted according to the procedure described in section Ul 
The full and dashed lines indicate a variation of the fragmentation parameter by ±lcr around 
the best fit value of a. The ratio R = MC/data is described in the caption of figure[5l 
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Figure 7: Normalised D*"^ meson cross sections as a function of zjct and Zhcm for the D*"^ jet 
sample. The data are corrected for hadronisation effects (see section [5]). The measurements 
are normalised to unity in the displayed range of zjct and z^cm, respectively. The data are 
compared to NLO predictions of HVQDIS with the Kartvelishvili parametrisation using the 
fragmentation function parameter a fitted according to the procedure described in section Ul 
The full and dashed lines indicate a variation of the fragmentation parameter by ±lcr around 
the best fit value of a. The ratio R = MC/data is described in the caption of figure[5l 
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Figure 8: Normalised D*^ meson cross sections as a function of Zhcm and zjct for the D*^ 
jet sample in two regions of W , for 14^ < 170 GeV and 11/ > 170 GeV. The measurements 
are normalised to unity in the displayed range of zjct and Zhcmj respectively. In addition, the 
MC predictions of RAPGAP are shown using the ALEPH setting for PYTHIA and the fitted 
fragmentation parameters a. for the Kartvelishvili parametrisation as given in table [8] for zjct 
and Zhcm- The ratios R = MC/data, as described in the caption of figure [5l are shown for both 
regions of W on top of each other. 
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Figure 9: Normalised D*"^ meson cross sections as a function of Zhom for the no D*"^ jet 
sample. The measurements are normalised to unity in the displayed range of Zhcm- The data 
are compared to MC predictions of RAPGAP with the ALEPH setting for PYTHIA and the 
Kartvelishvili fragmentation function using the fragmentation parameter a fitted according to 
the procedure described in section |7J The full and dashed lines indicate a variation of the 
fragmentation parameter by ±lcr around the best fit value of a. The dotted line shows the 
prediction of RAPGAP with the fragmentation parameter a = 4.4 extracted from the D*^ jet 
sample. The ratio R = MC/data is described in the caption of figure [51 
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Figure 10: Normalised D*"^ meson cross sections as a function of Zhcm for the no D*^ jet 
sample. The data are corrected for hadronisation effects (see section O. The measurements are 
normalised to unity in the displayed range of Zhcm- The data are compared to NLO predictions 
of HVQDIS with the Kartvelishvili fragmentation function using the fragmentation parameter 
a fitted according to the procedure described in section U\ The full and dashed lines indicate 
a variation of the fragmentation parameter by ±lcr around the best fit value of a. The ratio 
R = MC/data is described in the caption of figure[5l 
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